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Abstract 

A one dimensional kinetic Ising model at a finite temperature on a semi- 
infinite lattice with time varying boundary spins is considered. Exact 
expressions for the expectation values of the spin at each site are ob- 
tained, in terms of the time dependent boundary condition and the ini- 
tial conditions. The solution consists of a transient part which is due to 
the initial condition, and a part driven by the boundary. The latter is 
an evanescent wave when the boundary spin is oscillating harmonically. 
Low- and high-frequency limits are investigated with greater detail. The 
total magnetization of the lattice is also obtained. It is seen that for any 
arbitrary rapidly varying boundary conditions, this total magnetization is 
equal to the the boundary spin itself, plus essentially the time integral of 
the boundary spin. A nonuniform model is also investigated. 
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1 Introduction 



Dynamical spin systems have played a central role in non-equilibrium statistical 
models. The Ising model is widely studied in statistical mechanics, as it is 
simple and allows one to understand many features of phase transitions. The 
non-equilibrium properties of the Ising model follow from the spin dynamics. In 
his article 1 , Glauber introduced a dynamical model formulating the dynamics 
of spins, based on the rates coming from a detailed balance analysis. It is a 
simple non-equilibrium model of interacting spins with spin flip dynamics. An 
extension of the kinetic Ising model with nonuniform coupling constants on a 
one-dimensional lattice was introduced in [2] . In [3] , a damage spreading method 
was used to study the sensitivity of the time evolution of a kinetic Ising model 
with Glauber dynamics against the initial conditions. The full time dependence 
of the space-dependent magnetization and of the equal time spin-spin correlation 
functions were studied in [4] . Non-equilibrium two-time correlation and response 
functions for the ferromagnetic Ising chain with Glauber dynamics have been 
studied in [5, 6 . The dynamics of a left-right asymmetric Ising chain has been 
studied in [TJ. The response function to an infinitesimal magnetic held for 
the Ising-Glauber model with arbitrary exchange couplings was addressed in 
[8 . In [9], a Glauber model on a one-dimensional lattice with boundaries was 
studied, for both ferromagnetic and anti-ferromagnetic couplings. The large- 
time behavior of the one-point function was studied. It was shown that the 
system exhibits a dynamical phase transition, which is controlled by the rate of 
spin flip at the boundaries. 

It was shown in [TUl E] that for a nonuniform extension of the kinetic Ising 
model, there are cases where the system exhibits static and dynamical phase 
transitions. Using a transfer matrix method, it was shown that there are cases 
where the system exhibits a static phase transition, which is a change of behav- 
ior of the static profile of the expectation values of the spins near end points 
[10] - Using the same method, it was shown in that a dynamic phase tran- 
sition could occur as well: there is a fast phase where the relaxation time is 
independent of the reaction rates at the boundaries, and a slow phase where the 
relaxation time does depend on the reaction rates at the boundaries. 

Most of the studies on reaction diffusion models have been on cases where 
the boundary conditions are constant in time. Among the few models with time 
dependent boundary conditions, is the asymmetric simple exclusion process on a 
semi-infinite chain coupled at the end to a reservoir with a particle density that 
changes periodically in time 12 j. The situation is similar regarding the case 
of the kinetic Ising model as well. Among the exceptions are the study of the 
dynamical response of a two-dimensional Ising model subject to a square-wave 
external field |13) . and the study of a harmonic oscillator linearly coupled with 
a linear chain of Ising spins [21 [15] . 

In this article a one dimensional Ising model at temperature T on a semi- 
infinite lattice with time varying boundary spin is investigated. The paper is 
organized as follows. In section 2 a brief review of the formalism is presented, 
mainly to introduce the notation. In section 3, a semi- infinite lattice with 
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oscillating boundary spin is studied. The exact solution for the expectation 
values of the spin at any site is obtained. It is shown that there the boundary 
produces an evanescent wave in the lattice. The low and high frequency limits 
are studied in greater detail. The total magnetization of the lattice, M(t), is 
also obtained. It is shown that for rapidly changing boundary conditions, the 
total magnetization is equal to the the boundary spin itself, plus something 
proportional to the time integral of the boundary spin. A nonuniform model in 
also investigated. It is shown that its evolution operator eigenvalues are real. 
For the specific case of a two-part lattice with each part being homogeneous, 
the reflection and transmission coefficients corresponding to a harmonic source 
at the end of the lattice are calculated. Finally, section 4 is devoted to the 
concluding remarks. 

2 One-dimensional Ising model with nonuniform 
coupling constants 

Consider an Ising model on a one-dimensional lattice with L sites, labeled from 
1 to L. At each site of the lattice there is a spin interacting with its nearest 
neighboring sites according to the Ising Hamiltonian. At the boundaries there 
are fixed magnetic fields. Denoting the spin at the site j by Sj, and the magnetic 
field at the sites 1 and L by *8i and <8_l, one has for the Ising Hamiltonian 

L— 

H = - ^ s a-^ a+)I - ( BiSi- , Bl S L . (1) 

where J a is the coupling constant in the link a, and 

"4 (2) 

The link a links the sites a — \i and a + fi, so that a ± /j, are integers, and a 
runs from /i up to (L — fi). Throughout this paper, sites are denoted by Latin 
letters which represent integers, while links are denoted by Greek letters which 
represent integers plus one half (/j,). The spin variable Sj takes the values +1 
for spin up (j), or —1 for spin down (|). Define 

'/3J a , l<a<L 
K a :=\p<B u a = f i (3) 
J3K Ll a = L + [i 

where 

and fee is the Boltzmann's constant, and T is the temperature. Denoting the 
reaction rate from the configuration A to the configuration B by uj(A — > B), 
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and assuming that in each step only one spin flips, detailed balance demands 
the following for the reaction rates. 

u>[(S', Sj) -> (S", -Sj)} = Tj [1 - sj tanh(^_ M + K j+ „ s j+1 )}, 

Kj<L, (5) 

u[(S', si) -> (5', = Ti [1 - Sl tanh(^ + s 2 )], (6) 

w[(S ,/ , s L ) -> (5', -s L )] = Tl [1 - sl tanh(A L - M sl-i + Kl+„)]. (7) 

Tj's are independent of the configurations. For simplicity, we take them to be 
independent of the site. Then, rescaling the time they are set equal to one. 
So the evolution equation for the expectation value of the spin in the site j 

is 

— (sj) = -2 ( Sj ) + [tanh(^_ M + K j+fi ) + tanh(^-_ At - K j+Il )\ (sj-i) 

+ [tanh(-Kj_ M + K j+il ) - tanh(A' J _ At - K j+t ,)} (sj+i), 1 < j < L, 

^(si) = -2 (si) + [tanh(A M + K 1+fi ) + tanh(A M - K 1+Il )] 
+ [tanh(A M + K 1+IX ) - tanh(A^ - Ki +M )] (s 2 ), 

^(s L ) = -2 (s L ) + [tanh^L-M + K L+t ,) + tanh(A L - M - K L+li )] (s L -i) 

+ [tanh(Jf l-m + - tanh( A L - M - A L+M )] . (8) 

These can be written in the form 

^r{sj) = -2 (sj) + [tanh(iT 3 -_ M + if,- +|t ) + tanh(Aj_p - A^ +/ J] (s^-i) 

+ [tanh^-p + - tanh(^_ M - tf i+/t )] (s J+1 ), 1 < j < £, 

(9) 

(so) = 1, (10) 

(SL+l) = I- (11) 

3 Time varying boundary conditions on a semi- 
infinite lattice 

Consider a lattice for which the boundary spins (so and sl+i) are externally 
controlled, but the reactions at the internal sites satisfy detailed balance. The 
evolution equation is then the same as ((9]), but combined with boundary condi- 
tions different from (fit)]) and (TTTT) . A semi-infinite lattice the boundary of which 
is externally controlled, is obtained by letting L tend to infinity, and using the 
following boundary conditions 

<*o> - f(t), (12) 
(sj) does not blow up as j tends to infinity, (13) 
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instead of JTDJ and (|TT|) . 

A general solution of (|9]), combined with (|T2l) and (fl3l) . can be written as 
the sum of a particular solution plus a general solution of ©, combined with 
the homogeneous boundary conditions. 



3.1 Semi-infinite lattice with uniform couplings: the ho- 
mogeneous solution 

For a lattice with uniform couplings, K^s are denoted by K. The solution to 
the homogenous equation (vanishing /) is denoted by (fij)h- One arrives at 

— {s j ) h = -2{s j ) h +{tsnah{2K)]({s j _ 1 ) h + {s j+1 ) h ), < j. (14) 

Defining 

(Sj)h ■= -(s-j)h, J < 0, (15) 

one arrives at 

-< Sj ) h = -2 ( Sj ) h + [tanh(2 K)] « Si _i) h + (s j+1 ) h ), (16) 

which holds for all integers j. Denoting the linear operator acting on (s;)'s in 
the right-hand side of (|16l) by h, the above equation is of the form 

±( aj ) h = h) { Sl ) h , (17) 

where hj$ are the matrix elements of h. Defining the generating function G 
through 

oo 

G(z,t):= J2 zHsMt), (18) 

j=-oo 

one arrives at 

— = [-2+(z + z- 1 )tanh(2K)}G, (19) 

resulting in 

G{z,t) =exp{[-2 + (z + z- 1 ) tanh(2iT)]t}G(z,0), 

OO 

= exp(-2i) z k I k [2t tanh{2K)]G{z,0), 

k— — OO 

OO oo 

= exp(-2t) 7j 12 Ii-l[2ttanh(2A-)]<fl,) h (0), (20) 

j — — oo / — — oo 
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where is the modified Bessel function of first kind of order k. (f^O)) results in 

oo 

(sj) h (t) =exp(-2t) Ii-i[2*tanh(21O]<*i>h(0), 

/ — — OO 
OO 

= exp(-2t) ^2{lj-i[2t tanh(2 #)] - l j+ i[2t tanh(2#)]} (s/) h (0). 

(21) 

Using the large argument behavior of the modified Bessel functions, it is seen 
that 

(sj)h(t) ~ exp{-2 [1 - tanh(2 K)] t}, j > 0, (22) 
showing that the homogeneous solution tends to zero at large times. 

3.2 Semi-infinite lattice with uniform couplings: the par- 
ticular solution corresponding to harmonic boundary 
conditions 

The harmonic boundary condition is 

(s ) = Re[cr exp(— iuit)] (23) 
The following ansatz for a particular solution (sj) p to equations (J9)) and (|T2|) 

(sj)p = Re[crj exp(— iwt)] (24) 

results in 

(iw- 2)<rj + [tanh(2 J fC)](cr i+ i + cxj-i) = 0. (25) 
This has a solution of the form 



Tj = cz J , 



where z satisfies 



z + z 1 = 



(26) 



(27) 



tanh(2.?0 ' 

It is obvious that changing the sign of K results in changing the sign of z, while 
changing the sign of u> results in changing z to its complex conjugate. So it 
is sufficient to consider only nonegative values of K and u. From now on, it 
is assumed that K and co are nonnegative. (|27j) has two solution for z, which 
are inverse of each other, and none are unimodular. The boundary condition 
at infinity imposes that of the two solutions of type (|26p , only that solution is 
acceptable which corresponds to the root of (|27|) with modulus less than one. 
From now, only this root is denoted by z: 

z := r cxp(i6>) (28) 
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where r and 9 are real and r is positive and less that one. The solution to (|25p 
is then 

o =oqzK (29) 

As \z\ is less than one, the particular solution (IM|) describes an evanescent wave. 
Obviously, the rate of decay length and the phase speed, I and v respectively, 
satisfy 

mr 

v=-. (30) 

As the homogenous solution (|2"Tj) tends to zero for large times, the particular 
solution (|24[) is in fact the large times solution to the problem of harmonic 
boundary condition. 
Defining 



UJ 



tanh 2 K, 

r + r _1 



2 

the real and imaginary parts of (f27|) read 



(31) 



u cos 9 = — , 
b 

vV - 1 sin9 = |. (32) 

So u satisfies 

6 2 u 4 - (a 2 + b 2 + 1) u 2 + 1 = 0, (33) 
from which one arrives, for the solution which is larger than one, at 

1/2 

! ■ (34) 



1 + a 2 + b 2 + + a 2 + b 2 ) 2 - Ab 2 



\ 262 / ' 

This is increasing with respect to a, and decreasing with respect to b. Noting 
that 

du 1 ( 1 \ 

1 - -2 ) » (35) 

which shows that u is decreasing with respect to r, it is seen that r is decreasing 
with respect to w, and increasing with respect to K. One also has 

r = u — y/vF~—l. (36) 
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Regarding 6, differentiating the first equation in (j3"2")) with respect to a, one 



has 



resulting in 



du 89 
cos f — u sin — = 0, 



da 

dO cos 9 du 
da u sin 9 da ' 



9a 



(37) 



suit 



buy/(l + a 2 + 6 2 ) 2 -46 2 a 
fe 2 u 2 — 6 2 sin 6> 



Equation (|34|) shows that 
from which it is seen that 



bu(b 2 u 2 — u 2 ) a 



bu > 1, 



„ d9 sin0 

0< 7T ^ ■ 

aa a 



(38) 



(39) 



(40) 



The first inequality shows 9 is an increasing function of a, so it is an increasing 
function of oj. The second inequality results in 



d9_ < 9_ 

da ~ a : 



(41) 



which shows that (0/a) is a decreasing function of a. So (a/ 9) is an increasing 
function of a, or (uj/0) is an increasing function of uj. 
One also has 



d(2 b 2 u ) 

a^ 2 



= 1 



^(l + a 2 + 6 2 ) 2 -4& 2 



(42) 



and as 



y/{l + a 2 + b 2 ) 2 -4b 2 > 1 - b 2 , (43) 

it turns out that (bu) is increasing with b, so that is increasing with b. Hence 
(oj/9) is decreasing with K . 

The asymptotic behavior of r and 9 is summarized as 



tannic, 
tanh(2 K) 

5 

tanh(2if) 

VITTJ 2 ' 



V8 + w 2 + w Vl6 + cj 2 - V^ 2 + cj\/16 + cj 2 

71 : 



uj < 1, 
1 < cj, 

if« 1, 
1 < if 



(44) 
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tanh(2 K) 



Figure 1: 

The plot of r versus to for different values of tanh(2 K) 



and 



to cosh(2 K) 
2 • 

7T 

2' 

tan — , 



cos 



' 8 + to 2 - to V16 + to 2 



w < 1, 



1 < -FT 



(45) 



Among other things, it is seen that the phase speed, at low frequencies ap- 
proaches the constant value 2/[cosh(2AT)], while at high frequencies varies like 
(2to/n). 

Figure Q] is a plot of r versus to for different values of t&nh(2 K) from 0.1 to 
0.9. Figure [5] is a plot of the phase speed (to/ 9) versus to for different values of 
tanh(2K) from 0.1 to 0.9. 

The total magnetization, defined as the sum of the expectation values of the 
spins, is denoted by M. At large times only the particular solution contributes 
to the magnetization. So, 



M = Re 



CO 



1-z 



exp(— itot) 



t — >• oo. 



For the time- independent boundary condition, this leads to 

00 



M 



1 - tanh K ' 



(t oo, w = 0). 



(46) 



(47) 
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(u/9) 



o.i 



tanh(2 K) 



0.9 




Figure 2: 

The plot of the phase speed (oj/0) versus lo for different values of tanh(2i(f) 



For high frequencies, 



M = Re I o-Q 



tanh(2 K) 



-1U 



exp(— \ujt) > , (t — > oo, w — > oo). (48) 



This can be simplified to 



M = Re <^ ctq 



1 + 



tanh(2 if) 



—ico 



cxp(— iujt) > , 



= (s )(t) + [tanh(2 JT)] [S (t) - S ], (t -> oo, w -> oo), 



(49) 



where 



So(*) := / dt , <*o>(t'). 
Jo 



Sq := lim 



f / dtS (t) 



(50) 



One then arrives at a similar result for the magnetization when the boundary 
condition is any arbitrary rapidly varying function of time (so that its low 
frequency components are negligible): 



M = (« )(t) + [tanh(2J0] [S (t) - S ], 

(t — »■ oo, rapidly varying boundary conditions). 



(51) 
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3.3 Semi-infinite lattice with two parts of uniform cou- 
plings: the particular solution corresponding to har- 
monic boundary conditions 

Consider a semi-infinite lattice consisting of two parts, so that 

K a = \ KU a<N > (52) 
[K 2 , a > N K ' 

The time evolution equations for the expectation values of the spins are 

(ij) = -2 ( Sj ) + [tanh(2if 1 )] + < j < N, (53) 

(s N ) = -2 (s N ) + k_ (sjv-i) + k+ (sn+i), (54) 

(sj) = -2 ( Sj ) + [tanh(2 tf 2 )] « Sj _x) + (s j+1 )), N < j. (55) 

where 

k_ := tanh(if x + K 2 ) + tanh(ivT 1 - K 2 ), 

k + := tanh(ifi + K 2 ) - tanh(i^i - K 2 ). (56) 

Applying a harmonic boundary condition (|2"3"|) . one has for the particular solu- 
tion of the kind (|2"4"|). 



{A lZ {+B x z^), 0<j<N 
A 2 zi, N<j 



(57) 



where 

Zl + z ^-^ky / = 1 ' 2 ^ 

and \zi\ is smaller than one. The boundary condition results in 

A±+Bi = a . (59) 
From (j57| for j = N, one arrives at 

A 1 z? + B 1 Zi N =A 2 z?. (60) 

Finally, (|54p results in 

K _ (At zf- 1 + Si zf w+1 ) + k+ v4 2 z^ +1 = (-iw + 2) A 2 zf. (61) 
Equations (f59l) through (jBTT) give 





(«- 


Zl - 


f K + Z2+iw-2)Z[ 


^0 






hi _ 




_1 ) + (k+ z 2 +ioj — 


2) (zr w - 


zf) 




"(«- 


- z r 


1 + k + z 2 + iuj - 2) 


zf a 






hi _ 


K. 


_1 ) + (/c+ z 2 +iw- 

- (21 - z r x ) °o 


2)(z^- 


zf) 



A x = 
Bi = 

2 «_ (zr w+1 - zf- 1 ) + (« + z 2 + i W - 2) ( Z { N -z?y 



ii 



For large N, these simplify to 
Ax = (T , 

— (k_ z^ 1 + k + z 2 + iuj — 2) N a 



Bi = 



K_ Zl + K+ Z2 + 1 W 



K_ (z-r — Zi 1 ) Z Q ^ Z^ (Tn 

A 2 = — 1 ' 2 . 1 ° 63 

K_ Zl + K + Z2 + 1W — 2 

It can be easily shown that for the nonuniform lattice and at high frequencies, 
up to first term in cu^ 1 the magnetization is similar to the case of the uniform 
lattice. 



3.4 Semi-infinite lattice with nonuniform couplings: the 
relaxation times 

The general solution of ([§]) with (fT2"j) and (fT3"|) is the sum of a particular solution 
and the general solution to © and (fT2"|) and (fTBf with vanishing /. The latter 
(the homogeneous solution) satisfies 

^■(sj)h = -2 (sj)h + [tanh(-Kj_ M + if,- +Al ) + tanh(i^_^ - (sj-i)h 

+ [tanh(ifj_ M + - tanh(if J _ fl - (sj+i)h, 1 < 3 < L, 

(*o)h = 0, 

(sl+iW = 0, (64) 

which can be written as (|17[) . Denoting an eigenvalue of h by £7, and the 
corresponding eigenvector by -0b, it is seen that there are solutions to of 
the form 

<*j>h(*) = 1>Ej exp(Et). (65) 
These solutions decay with a relaxation time r satisfying 

1 



Re(£) ' 



(66) 



One can see that the eigenvalues of the operator h are real. To see this, one 
notices that equations (|64|) are the same as the equations corresponding to the 
homogeneous solution of (JHJ)- So the homogeneous solution to the Ising chain 
externally driven at ends, is the same as the homogenous solution to the Ising 
chain with magnetic fields at boundaries. The evolution equation for the latter 
satisfies the detailed balance. For any evolution satisfying detailed balance, the 
eigenvalues of the evolution operator are real. To see this, one notices that the 
criterion of the detailed balance is 

u(B -> A) = Yg exp[/3(£ B - £ A )], (67) 

where A and B are two different state, Yg s are real nonnegative numbers (for 
B 7^ A) satisfying 

Yf = Y£, (68) 
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and £ is the energy of the system in the state A. So the matrix Y is Hcrmitian. 
Equation (|6T[) means that the evolution matrix H, the components of which are 
uj(B — > A)'s, is a similarity-transformed of Y. As Y is Hermitian, the eigenvalues 
of Y are real. As H is a similarity-transformed of Y, the eigenvalues of H are 
the same as the eigenvalues of Y . So the eigenvalues of H are real ([TB] for 
example). The eigenvalues of h are eigenvalues of H as well. So the eigenvalues 
of h are real. 

4 Concluding remarks 

A one dimensional kinetic Ising model at temperature T, with time varying 
boundary conditions was studied, for the case the lattice is semi-infinite. The 
evolution equation for the expectation values of the spins was investigated. For 
the case of harmonic boundary conditions, with uniform couplings, exact par- 
ticular solutions were obtained for the expectation values of the spins, as well 
as the total magnetization. The low- and high-frequency behaviors were stud- 
ied in more detail. Models for which the coupling constant is nonuniform were 
also studied. Physically, such nonuniform couplings could arise when either the 
interaction between spins or the temperature depends on the position. As a spe- 
cific example, the harmonic solution on a semi-infinite lattice consisting of two 
homogeneous parts studied. Finally, it was shown that for a general (nonuni- 
form) lattice, the eigenvalues corresponding to the evolution operator are real. 

Acknowledgement: This work was supported by the research council of the 
Alzahra University. 
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